Loss of control in pattern-directed nucleation: a theoretical study 
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The properties of template-directed nucleation are studied close to the transition where full nu- 
cleation control is lost and additional nucleation occurs beyond the pre-patterned regions. First, 
kinetic Monte Carlo simulations are performed to obtain information on a microscopic level. Here 
the experimentally relevant cases of ID stripe patterns and 2D square lattice symmetry are consid- 
ered. The nucleation properties in the transition region depend in a complex way on the parameters 
of the system, i.e. the flux, the surface diffusion constant, the geometric properties of the pattern 
and the desorption rate. Second, the properties of the stationary concentration field in the fully 
controlled case are studied to derive the remaining nucleation probability and thus to characterize 
the loss of nucleation control. Using the analytically accessible solution of a model system with 
purely radial symmetry, some of the observed properties can be rationalized. A detailed comparison 
to the Monte Carlo data is included. 
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I. INTRODUCTION 

The controlled fabrication of nano-structures on sur- 
faces has achieved a lot of interest during the past decade 
because of its great application potential in the area of 
micro-electronics. Predefined structures at this scale can 
in principle be obtained in a number of ways: by self- 
organizing growth, by photo-lithography and etching, 
and by template directed nucleation. The occurrence of 
self-organizing growth mechanisms is relatively rare and 
the obtained patterns are by no means arbitrarily control- 
lable; the process of lithography and etching has become 
very sophisticated and versatile, but unfortunately it is 
not applicable to every kind of substrate; many organic 
substances for example cannot be patterned in this way. 

This is the point where template directed growth be- 
comes interesting: in this technique the areas of desired 
accretion are marked in some way to favor the nucleation 
of the substance, which is subsequently deposited. This 
way, the patterning is achieved indirectly by means of an 
already existing technique. The preferential nucleation 
sites can be created in numerous ways: the flat surface 
can be roughened at the points of interest by techniques 
such as focused ion beam engraving 1 , laser ablation 2 , 
direct imprinting 3 , or even by the tip of atomic force 
microscopes^; alternatively substances with a compara- 
tively higher binding energy to the depositable substance 
can be applied in a controlled fashion to provide a pre- 
defined nucleation site. This can be accomplished e.g. 
by printing^, by dip-pen writing?, or by lithographic 
methods 8 . 

The specific process analyzed in this work is that of 
template-directed nucleation as applied by Wang et al£ 
However, due to the general nature of the underlying 
mechanisms, the conclusions drawn in the course of the 



paper should also be applicable to a wide range of similar 
techniques. In their experiments a silicon oxide surface 
was patterned with gold dots electron beam lithography. 
Then organic substances with a comparably higher affin- 
ity to the gold pattern were deposited by molecular beam. 
Following adsorption these molecules diffused along the 
surface to eventually either nucleate on a gold dot, to 
form a new two-dimensional nucleus in combination with 
other adatoms, or to desorb as visualized in ngureQ] The 
experiments were carried out at varied pattern spacing 
and for several organic substances, of which the deposi- 
tion of diferrocene exhibited the least fluctuations. 

From an application point of view one is interested in 
a maximal range of full nucleation control. Specifically 
this means: how far can (at constant deposition flux and 
temperature) the pre-patterned gold structures be placed 
apart until additional nucleation will occur aside of the 
pattern? This work focuses on the theoretical under- 
standing and description of this transition region, which, 
as a first approximation, one expects to occur, when the 
density of the pattern matches the saturated island den- 
sity that would be obtained from a nucleation experiment 
on an unpatterned surface under the same conditions 1 . 
On second thought, however, a sensible and regular ar- 
ray of sinks should provide a better 'drainage' of the de- 
posited material and hence an increased nucleation con- 
trol. 

We approach this problem by means of two comple- 
mentary theoretical methods. In section |TT] we present 
a detailed description of the MC simulation applied to 
model the experiment. Then, in section [TTTI introduce the 
dimensionless notation from Ref. @, which will be applied 
here, as well, and we discuss the qualitative behavior of 
the nucleation control. 

To obtain a theoretical understanding of the limit of 
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nearly full nucleation control, in section IIVI we apply a 
model similar to the level-set methods of Ref. l9l lldfil"l . 
In this approach the deposited adatoms are assumed to 
diffusive independently such that their concentration can 
be described in BCF-likc 12 fashion by a continuous con- 
centration field, to which the pre-patterned sites pro- 
vide the level-set boundaries. The nucleation probabil- 
ity is then derived on the basis of standard nucleation 
theor y 13 i 14 ' 15 i 16 . Specifically, the local nucleation rate is 
a function of the adatom concentration and hence the 
overall nucleation probability can be expressed by a cor- 
responding integral over the surface. On this basis, the 
nucleation behavior is then discussed in the one- and 
two-dimensional case analytically as well as numerically. 
As the deposition flux is experimentally easily accessible, 
this is done with special respect to the flux dependence. 



which the interaction energy is zero and interactions can 
be neglected: 




FIG. 1: Schematic representation of the considered mecha- 



II. MONTE CARLO SIMULATIONS 

To gain deeper insight on the atomistic mechanisms of 
pattern directed growth we examine the behavior of pat- 
terned systems by means of a Monte-Carlo (MC) simula- 
tion. Since we are not interested in the short-time ballis- 
tic movement of the adsorbed particles, a MC algorithm 
is sufficient to model the diffusive behavior-il. It has been 
show n 18 ' 19 that a MC step of a diffusive process can be 
associated to a corresponding mean time step, allowing 
a dynamical interpretation. In contrast to the alterna- 
tive of a molecular dynamics (MD) simulation, a MC 
algorithm provides the additional possibility of enhanced 
computational efficiency due to exclusion of 'spectator- 
particles' from the propagation routines, the substrate 
particles in this case. 

Since the molecular dimensions are much smaller than 
the relevant environmental length scales we use simple 
Lennard-Jones (LJ) pair potentials to mimic the molec- 
ular interactions. The parameters a and e are a measure 
for the atom-atom distance and depth of the potential 
well, respectively. The substrate (S), the pattern (A), 
and the deposited material (M) are taken into account 
as individual particles. To enhance computational effi- 
ciency, the chosen LJ variant possesses a cutoff r c , above 
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The additional terms ensure a continuous transition of 
potential and force at r c while showing only a minor in- 
fluence on the potential otherwise. 

With all of the A and S particles fixed, the relevant 
interactions are reduced to the ones involving M par- 
ticles. The chosen LJ parameters are ema/^bT = 5, 
emm/^bT = 3.7 and cms/^bT = 1.5, while for the sake 
of simplicity all radial constants are set to a = 0.890 a, 
with a being the fee lattice constant and r c = 2.5 a. The 
energy parameters were chosen such as to prefer the ac- 
cretion of M on A, as well of M on M over the simple 
adsorption on the substrate. The M-S interaction is a 
compromise between a sufficiently high surface diffusion 
rate on the one hand and at reasonably low desorption 
rate on the other. 

We apply off-lattice displacements. Per sweep every 
particle is subjected to a trial move, which is a product 
of a randomly generated vector and an also stochasti- 
cally determined jump length. The normalized displace- 
ment vector is chosen from a uniform distribution on a 
sphere^, while the jump length is a uniform deviate in 
the interval [0.01 a, 0.5 a]. On the cost of computational 
efficiency, this interval is chosen small enough to prevent 
frequent 'tunneling' of particles through energetic barri- 
ers. 

The simulation is conducted in a simulation box with 
periodic boundary conditions. A typical simulation setup 
consists of a substrate located at the middle of the z-axis 
covering the xy-plane, yielding a quasi- infinite surface. 
Since the chosen potential is radially isotropic the possi- 
ble crystal lattices are fee or hdp. However, neglecting 
the proper structure of the material should not pose any 
problem: firstly, because simple lattices have previously 
been successfully used to model the growth of different 
materials^; and secondly because all necessary physical 
properties - especially the preference of the M-A interac- 
tion over the M-S alternative - are included. The lattice 
of the substrate is oriented with (lll)||z. The 'unit- 
cells' resulting from this orientation have the dimensions 
of a x V3 a x v6 a and contain six particles in three layers. 
The substrate is then patterned with three more layers 
of A particles in the desired morphology. As above, the 
same lattice constant is used for both, S and A particles. 
A typical setup resembling the experimental conditions 
of Ref. d is depicted in figure O 

The deposition flux is modeled by periodic addition of 
M particles to the substrate. The x- and y-coordinate of 



3 




FIG. 2: Typical simulation setup for a pre-patterned surface 
with dots on a squared lattice. The substrate has the dimen- 
sions of 70 x 35 unit-cells. The core of the middle cluster is 
built from pattern material depicted in white. 




p/2 2 s p/2 



FIG. 3: A typical simulation setup consisting of a substrate 
with dimensions 50 x 20 as well as a pattern stripe on the 
substrate. Close to the middle between the stripes - with 
respect to the periodic boundary conditions - a new nucleus 
is formed. 



remain with the presented parameters. 

The desorption coefficient was determined from a simi- 
lar simulation: after equilibration adsorption and desorp- 
tion flux are equal requiring F — g c. The desorption co- 
efficient g = 12.7 • 10~ 6 /step could hence be derived from 
the average equilibrium adatom concentration. Conse- 
quently, the influence of desorption cannot be neglected 
with respect to the stationary concentration (for more 
details see equation ([25])). These parameters and a mean 
system size of p = 40 a in combination with the consider- 
ations of appendix [5] equation (|A4[) lead to an estimate 
of 50 000 steps to reach a stationary concentration profile 
in the ID-case. This lies well within the typical simula- 
tion time of 2 — 4 • 10 6 steps and is sufficiently below 3 • 10 5 
steps, the time at which the substrate becomes saturated 
with islands. 



III. NUCLEATION CONTROL 

We regard the loss of nucleation control an a quadratic 
pattern as experimentally observed in Ref . |8| with special 
respect to the transition region, where control is gradu- 
ally lost. To this purpose we simulated substrates con- 
structed of n x n/2 x 1 'unit-cells', where the ratio of x- to 
y-dimension is chosen to approximate the desired square 
lattice by means of the non-square 'unit-cells'. A dot of 
radius s — 3a and three layers of height built from the 
same lattice are added to the middle of the surface. 



the deposition point are determined uniformly within the 
dimensions of the simulation box. The particles are then 
inserted at maximal z and subsequently moved toward 
the substrate, until an energy minimum is met. After 
this, the particles are subject to the described MC prop- 
agations. Desorption is included as well : the system is 
frequently probed for a particle out of contact (r > r c ) 
with any other particle. If such a particle is found, it is 
considered desorbed and removed from the simulation. 

The simulations, with the exception of the latter ones 
of section IIV C 21 were conducted at a deposition rate 
of F = 1.43 ■ 10~ 7 pt/(<7 2 step), where pt = particles 
and 'step' refers to a MC-step. The diffusion coefficient 
was determined by analysis of the mean square displace- 
ment of non-interacting M particles (6mm = k^,T) on 
an unpatterned surface and found to be D — 3.50 • 
10- 3 er 2 / stc P- With a resulting rate of D/F = 2.45 • 
10 4 cr 4 /pt, which corresponds to D/F = 1.15 • 10 4 in the 
units of Brunei, the considered setup lies one order of 
magnitude below the condition of D/F > 10 5 assumed 
there for the theoretical treatment of non-interacting par- 
ticles. This in addition to the inclusion of desorption 
which is also typically neglecte d 13 ' 14 i 15 ' 16 i 21 can lead to 
deviations from the ideally expected behavior. As, how- 
ever, an increase of the temperature would not only raise 
D but desorption as well, and a sufficient decrease of 
F would significantly reduce the simulation efficiency we 



A. Dimensionless units 

To evaluate the obtained data and to allow the di- 
rect comparison of different experimental as well as sim- 
ulated systems, a dimensionless representation in analogy 
to Ref. [8| is introduced. We define the dimensionless nu- 
cleation efficiency as the number of prepatterned islands 
divided by the number of the actually obtained nuclei: 
for a single unit-cell of a square pattern the nucleation 
control is hence given by 



where R is a fractional number reflecting the averaged 
number of additional nuclei per unit cell. 

The spatial coordinate is renormalized by the charac- 
teristic length scale A, which is derived from the nucleus 
density of an unpatterned surface: 




where A is the surface area. Using A the lattice constant 
p of the pattern can be expressed in the dimensionless 
form: 
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It should be noted that for patterns with more than one 
lattice constant the specific choice of p is in principle 
arbitrary. 

B. Limit behavior 

Mechanistically, one can expect growth to occur by a 
diffusion-aggregation mechanism as depicted in figure [T] 
After deposition, the adatoms explore the surface diffu- 
sively, until they meet or form a nucleus or pattern site, 
or desorb from the surface altogether. 

On an unpatterned surface nucleation occurs if a suffi- 
cient number of particles aggregate to form a stable two- 
dimensional nucleus, that will then grow by accretion of 
further particles. On a pre-patterned surface on the other 
hand, the first nuclei are already supplied by the pattern 
material, which is chosen for its comparatively more fa- 
vorable interaction energy with the deposited substance. 
The subsequent growth process is then an interplay of de- 
pletion, diffusion and deposition/desorption. In the limit 
where p is sufficiently small or D/F accordingly large, all 
arriving particles are transported to the nuclei and full 
nucleation control is maintained^. One consequently 
expects 

dn = 1 for p <C 1. (5) 

Should however the adatom concentration at some 
point on the surface facilitate the formation of additional 
nuclei control is gradually lost. This leads to the second 
limit case: at very large spacing between the dots the ex- 
istence of the pattern is negligible and the obtained nu- 
cleus density approaches the density of an unpatterned 
system. As the unpatterned system is used for the renor- 
malization of the space coordinate, the expected limit 
behavior is 

x N = p~ 2 for p > 1. (6) 

C. Transition region 

As a first approximation, one expects the loss of nucle- 
ation control to set in, when the density of the patterned 
surface matches that of its unpatterned equivalent. In 
the case of the applied quadratic pattern this corresponds 
to p = A and consequently p — 1. A sensible array of 
sinks should however supply a more effective 'drainage' 
of the deposited material, resulting in prolonged nucle- 
ation control for values of p > 1 . 

To further investigate the behavior at the transition 
point on the basis of the results in Ref. [H, we conducted 
additional simulations at different p. The first four sys- 
tems following the loss of nucleation control were ob- 
tained from twelve independent runs while the remaining 
points were determined from six runs. The number of is- 
lands is determined from the stationary regime before 
the onset of coalescence. The obtained data is depicted 



in figure [4] as circles exhibiting an extrapolated loss of 
nucleation at p = 1.27 with A = 24.8 a. 

It can be seen that the first three data points beyond 
the loss of nucleation control follow in almost linear suc- 
cession, whereas a jump is observed with respect to the 
following points, which are closer to the expected limit 
behavior. The standard deviation of these results, which 
is about the size of the symbols, suggests that this is not 
a statistical effect. Closer scrutiny of this jump reveals 
that it is observed as drops below 1/2, marking the 
point where a second additional nucleus is formed and 
which is indicated by the dashed line. For in > 1/2 only 
one additional island is formed which is expected to be 
found close to the corner, the spot of highest concentra- 
tion and farthest away from the pattern. Any additional 
nucleation will therefore in principle lead to a new lattice 
constructed from squares with an edge length smaller by 
a factor of 1/ y/2. This way the area is shared by the is- 
lands still quite 'effectively'. For % < 1/2 however this 
behavior breaks down, because the nucleation of the sec- 
ond island discontinues this scheme and approaches 
the disordered limit behavior. This becomes more so as 
nucleation control is lost further. 
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FIG. 4: Nucleation control in dependence of the reduced 
periodicity for the fluxes 1.43 • 10~ r pt/(step a 2 ) (o) and 
1.43 ■ 10 -6 pt/(step a 2 ) (x). As before the limit behavior for 
p < 1 and p S> 1 is indicated by dotted lines. The dashed- 
dotted line marks £n = 1/2. The dashed lines symbolize 
the deviation from the limit behavior caused by the dot size 
according to equation (Q) with s e g = 0.23 and 0.43. 

As with growing p the nucleation besides the pattern 
site becomes more and more statistical, the nucleation 
control can be approximated by the follwing considera- 
tion: taking a dot to occupy an area of its 2 in addition to 
its mean capture area of A 2 , and assuming the remaining 
area to be filled in the disordered nucleation limit one 
obtains 

«N = (p 2 - 7T S 2 ) 1 , (7) 

where the reduced radius is defined as s := s/A. However, 
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as s is only of the order of a few particles, the initial 
accretion of molecules around its boundary introduces 
a non-negligible size increase. Hence one has to apply 
an effective radius which can be derived from figure [7] of 
section HV C 21 one obtains s c ff = 0.23 and 0.43 which in 
combination with (O leads to the dashed lines in figure |U 
It can be seen that this expression indeed describes the 
convergence to the limit behavior, which occurs slower 
with decreasing A. 

From the experimental point of view, however, one is 
mainly interested in two questions: (a) how does the 
point at which nucleation control is lost (p) depend on 
the deposition flux (F)? And (b) how does the mean 
number of additional nuclei (R) change with increasing 
distance between the predefined dots (p)? These ques- 
tions are approached from a theoretical point of view in 
the next section. 



IV. THEORETICAL APPROACH OF THE 
TRANSITION REGION 

For the description of the transition region we adopt 
the following probabilistic picture. We consider a pattern 
which at given p, s, F, D and g is under (nearly) full nu- 
cleation control, i.e. « 1. In addition, we neglect the 
spatial growth of the pattern due to particle accretion 
and we require the adatom concentration c(x) to reach a 
stationary steady state based on the following consider- 
ation: according to standard nucleation theor y 13 i 14 ' 15 i 16 
an unpatterned surface will (after it passed through a 
regime of rapid nucleation) saturate with respect to the 
nucleus density (diV/di 0); hence c(x) becomes quasi- 
stationary as well. Because of its good experimental re- 
producibility and near independence on time we select 
this regime as reference state. If one now considers a pat- 
terned surface of comparable nucleus density, this surface 
can also be expected to exhibit a stationary adatom con- 
centration. Some basic considerations on the time neces- 
sary to reach this state can be found in appendix [A"l 

However, although the stationary system is supposed 
to be under full nucleation control, there remains a finite 
chance for nucleation due to fluctuations in the adatom 
concentration. Specifically this means that in the con- 
sidered scenario any additional nucleation is not caused 
by the dynamic increase of the concentration field which 
eventually supercedes a threshold, but only by the re- 
maining small but finite nucleation probability of the 
static concentration field. This probability can be evalu- 
ated by methods of standard nucleation theory, accord- 
ing to which the local rate of formation of new nuclei 
depends (after application of the Walton relation^) on 
the concentration as 

^occ(xf+\ (8) 

where all clusters of size i* or less are considered ther- 
modynamically unstable while larger clusters are stable. 



In analogy to the level set methods of e.g. Refs. [9lfl0l we 
determine the rate of nucleation in one unit-cell from the 
corresponding spatial integral. The mean number of ad- 
ditional nuclei formed up to time t can hence be expressed 
by temporal integration. Neglecting the time necessary 
to reach the steady concentration profile we find the di- 
rect proportionality 

Roct - [ d^c(x) r+1 oc ^ • / dAc(x)'* +1 (9) 
Ju F Ju 

with U as the unit-cell and a proportionality constant 
independent of p, s, F, D, g and t. The latter part holds 
for states of identical coverage 9. 

As the limit of full nucleation control with x-^ = 1 (or 
R = 0) is not well defined in this probabilistic approach, 
we consider a critical R with < R <C 1 , which serves to 
describe the point of transition. One can now examine 
the dependency of p, s, F, D 7 and g at constant R: if for 
example the deposition flux F is increased slightly, how 
will the corresponding p have to change? It should be 
noted that in the limit i* — * oo only the point of max- 
imal concentration contributes significantly to R, which 
in turn leads to the same result as a deterministic crit- 
ical concentration approach, where new nuclei are con- 
sidered to be formed when a concentration threshold is 
superceded. 

The qualitative flux dependency is intuitively accessi- 
ble: An increase in F results in a decrease of p and vice 
versa, but the behavior of p also depends on the density 
of the unpatterned substrate via A. Thus, if p and A ex- 
hibit different dependencies on F/D, or if the behavior of 
the pattern is influenced strongly by the dot diameter, it 
is conceivable that the specific point at which nucleation 
control is lost can be shifted. This behavior is examined 
in the following section. 

A. Nucleus densities 

To predict the point at which nucleation control breaks 
down we have to find expressions for the nucleus density 
on unpatterned and patterned surfaces as a function of 
the growth parameters. 

For two-dimensional nucleation on an unpatterned sub- 
strate standard nucleation theory-1^ supplies a general 
model, which is based on rate equations and leads to 
one of the central relations, namely 

fD\ x/2 

Aoc f -J with 1/3<x<1 (10) 

at constant coverage 0. This in turn supplies the de- 
sired F and D dependence of A, while the exponent is a 
function of the critical nucleus size i* according to 

x = -^—- (11) 

A i* + 2 y ' 

In the next step we investigate the nucleation behav- 
ior of patterned surfaces. Similar to BCF-theoryi^ we 
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apply the following continuum approach: neglecting the 
stabilizing effects of unstable clusters we seek the station- 
ary adatom surface concentration satisfying the diffusion 
equation 



flc(x) 
dt 



= DV 2 c(x) - g c(x) + F = . 



(12) 



On the basis of this stationary concentration field one can 
then obtain the mean number of additional nuclei accord- 
ing to equation ©. The dots are modeled by localized 
ideal sinks supplying c| s j n k = as boundary conditions. 

The general shape of the concentration field can be 
discussed by means of a reduced representation. Intro- 
ducing dimcnsionlcss scales for length and concentration 



p z b 



as well as the dimensionless parameter 



13:= ^fg~/Dp, 
equation (I12|) can be rewritten as 
V 2 c(x) - /3 2 c(x) = 



(13) 



(14) 



(15) 



supplemented by the corresponding boundary conditions. 
One can now easily see that any solution fulfilling equa- 
tion (fT2|) has the general form of 



F I3 2 

c(x) = - — c(x,/3). 
9 2 



(16) 



Thus, for g ^ the general shape of the concentration 
field, i.e. c(x), is completely determined by f3(g,D,p) 
in combination with the boundary conditions, while for 
g = it depends only on the boundaries. The absolute 
concentration, however, is always directly proportional 
to F. 

We now turn to the desired relation between p and F: 
demanding constant R as given in ^ and substituting 
the scaling behavior of (p~6j) one obtains 



RocF 1 " ■p 2 ^+ 1 )+ d - / dic(x,/3) 



(17) 



where d = 1,2 is the dimensionality of the considered 
system. Since in general the above integral has a complex 
dependence on p via f3, equation (|17[) cannot be written in 
the form p = h(F). Nevertheless important information 
can be gained from the logarithmic derivative in analogy 
to equation (fTU|) . We define 



9 log (FY 
9 'log (p) / 



1 



2 a* 



i) 



d 



log / dic(x,/3) l * +1 



(18) 



(19) 



with the help of which the behavior of the transition 
point can be expressed as 



p ex F^ 2+l ' a . 



(20) 



As the important information is contained in the expo- 
nent we additionally introduce 



X 
2 



(21) 



which is discussed in detail by means of numerical as well 
as analytical methods in section HV C II 

Analogously, the relation between R and p at constant 
flux is obtained from 



a log (R) 

d log (p) J p 
i* a . 



(22a) 
(22b) 



This relation provides a useful means of comparison to 
the simulation data discussed in section [TV C 21 with spe- 
cial respect to figure [H 



B. One-dimensional pattern 

1. Analytical results 

The one-dimensional case can be interpreted as a pat- 
tern of parallel stripes with spacing p between them and 
a width of 2 s as depicted in figure[3] In the dimensionless 
notation of equation (| 1 3j) the value of s is of no relevance, 
however, as the boundary conditions are supplied by the 
sinks with c(x = 0) = c(x = 1) = 0. For g = the sta- 
ble concentration profile is given by the simple parabolic 
equation 



;1D 



(x) 



which possesses a maximum of 



„1D 



p 2 F 
8D 



(23) 



(24) 



at x = p/2. 

Including desorption with g > identical boundary 
conditions lead to the stationary solution 



-ID 



= (- 7 exp [/?£]-(!- 7 ) exp[-/3£] + l) , (25) 



with (3 as defined in ((TU) and 7 := (1 + cxp[/3]) _1 . The 
maximum is found to be 



, , F / expf/3/21 
c(x = 1/2 = - 1 - 2 - V 'l 
9 V l + exp /3 



(26) 



In the limit g — > equations ([23)1 and (|2"6"|) correspond 
to (|2"3"|) and respectively, as can easily be seen from 
a Taylor expansion at p/2. With growing (3 higher order 
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terms begin to contribute, but below (3 ~ 1 their influ- i — i 1 1 1 1 • 1 1 r 

ence is rather small and the parabolic shape is retained 

in general. At further increase of (3 the function starts 0.25 - ( a ) P — 30 a ^ ^ 

to develop a plateau between the sinks. The maximal ' s s 

achievable concentration in the case /3 — ► oo corresponds 0.2 - / . ■ tj ■ "□' ■ □ n s 

to an equilibrium in the absence of sinks given by F/g or - / / ^ 

correspondingly 2I0 2 in dimcnsionlcss scales. ~oi5- ' • \ 

/ ••' ' n - \ 

Here, we refrain from an evaluation of equation (|18|) ' \ X N 

for the ID case, which can be easily done analytically 0-1 a \ 

for (3 = or numerically for (3 ^ 0. Rather, to get /•' '\ 

a better understanding of the stationary approximation, 05 ~ J? 

we analyze the concentration fields of the simulation. " r \ 

0.25 0.5 0.75 1 

X 



2. Simulation 



We simulated surfaces with a striped pattern to evalu- 
ate the validity of the mathematical considerations from 
section IIV Al The simulated systems consist of a sub- 
strate built from (30, 40, 50) x 20 x 1 unit-cells. The pat- 
tern stripe is constructed from the same lattice; it covers 
the y-dimension and possesses a width of 6 a along x. 
As can be seen in figure EJ the M particles preferably 
aggregate at or on top of the step. However, if p 2 F/D 
is sufficiently high, additional nucleation, usually close 
to the middle between the stripes, can be observed as 
expected. 

The mean adatom concentration along the x-axis was 
determined by an average over the full length of the sim- 
ulation time and y-axis. The results are given in fig- 
ure [3 At small distances between the stripes (p = 30 a) 
the obtained concentration dependence exhibits the pre- 
dicted behavior of non-interacting adatoms as given by 
equation ([25)1 (dotted line). When the steps are placed 
farther apart two effects can be observed: (a) the influ- 
ence of desorption on the concentration profile becomes 
more pronounced, which can be seen in comparison to 
the desorption-less concentration profile from equation 
(|23| (dashed lines); and (b) mutual stabilization of the 
adatoms becomes more important. At p = 40 a localized 
regions of higher concentration due to the formation of 
unstable clusters appear. This stands in good agreement 
with the visual observations frequently showing clusters 
of two or sometimes three atoms in the vicinity of the 
point of highest concentration. Finally, at p = 50 a nucle- 
ation sets in leading to a steep increase in concentration 
where c 3> 0.25. At this point it should again be noted 
that the averaged concentration profiles do not represent 
a single configuration but represent the mean of the sys- 
tem before and after nucleation. Additionally it should 
be kept in mind that only the reduced concentrations 
decrease with increasing p; the absolute values behave 
oppositely. 
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FIG. 5: Mean reduced concentration for different distances 
between pattern stripes. The dashed line represents the ideal 
behavior without desorption according to equation (|23[l . The 
dotted lines include desorption according to equation ((25}. 
For p = 30 a the results are close to the expected behavior. 
At p = 40 a stabilization of adatoms causes a domain of higher 
concentration close to the middle. At p = 50 a nucleation sets 
in resulting in a peak around x — 0.5 where c ^> 0.25. 



C. Two-dimensional pattern 

1. Analytical and numerical results 

In the one-dimensional case the concentration at the 
sink has to be zero and consequently the step width s 
does not have any influence on the concentration fields. 
In two dimensions however, the relation of sink size to 
pattern spacing is of importance. With the radial coor- 
dinate r and dot radius s as system parameters we again 
use dimensionless scales: 



f := r/p and s := s/p 



(27) 



Unfortunately, the analytic evaluation of c(x) is only 
available for a few limit cases. We therefore first adopt 
the approximation of a radially isotropic concentration 
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field, leading to the boundary conditions 

dc(f) 



c(r = s) = and 



df 



= 0. 



(28) 



In polar coordinates and without desorption equa- 
tion (fT2|) becomes 



1 d_ / dc 

r dr \ dr 



F 
D 



which after twofold integration yields 

l/.o .o\ 1 



c is °(f,i) = i(,s 2 -f 2 ) + 
with the maximum 



4 \ s 



CLc= i(-ln(2S) + 2S 2 -i 



(29) 



(30) 



(31) 



at f = 1/2. 

On the basis of this ansatz one can discuss the mathe- 
matically simple, but experimentally impractical desorp- 
tionless scenario where s and p are scaled varied propor- 
tionally, such that a = const: in this case the general 
shape of the stationary reduced concentration field re- 
mains unchanged as can be seen from equation (|16[) ; thus 
in combination with l|19p and (|2ip one obtains 



and 6 = 0. 



(32) 



The experimentally more realistic case where p is 
changed at constant absolute dot size s is analytically 
accessible in the following limits: first, in case i* — > oo 
only the point of maximal concentration contributes to 
R. One can thus apply the c max ansatz using the no- 
desorption approximation equation (|3ip . As shown in ap- 
pendix [B] this expression can be fitted to successfully de- 
scribe the numerically observed maximal concentrations. 
Combination of this expression with © and subsequent 
differentiation yields 



b = 



-1 + 4s 2 



2 (In (2) + In (a)) ' 



(33) 



This expression is included into figure U[a) as a thin dot- 
ted line. For s — > one finds b — > 0; as s increases 
the deviations from the assumed radial isotropy become 
dominant and the approximation loses its validity. 

The second analytically accessible setup is found in the 
limit (3 — > oo: In this case the influence of the sinks is 
diminishingly small and every part of the surface that is 
covered by pattern material has an adatom concentration 
of F/g. The integrated nucleation probability is then 
given by 



R oH [ — ) (p 2 — 7T , 

\9 ' 



(34) 



and differentiation according to (fT8)) and subsequent sub- 
stitution (|2"Tj) leads to 



b=\ 
2 vi* 



1 



(35) 



This relation is depicted in figure[6jd) for the correspond- 
ing i* by the thin dotted lines. 

Alternatively the stationary concentration for arbi- 
trary setups of lattice and dot shape is readily acces- 
sible by means of numerical solution: expanding the con- 
centration field to a discrete grid and expressing equa- 
tion (fT2|) by explicit finite differences, one obtains a sys- 
tem of linear equations of the form 



(D - g 1) c s 



(36) 



where D and f describe the concentration change with 
the next time interval due to diffusion and deposition 
flux, respectively. This set of linear equations can then 
be solved by standard matrix algebra 2 ™ On the basis 
of these solutions one can now determine the scaling be- 
havior of the nucleation in the following way: first, the 
patterned exponent a is determined by numerical differ- 
entiation according to (|18[) at constant s, D and g by 
variation of p by the width of one grid point. However, 
the overall curves are calculated at constant (3: when at 
constant s the distance p is increased to realize a different 
s, g is adapted appropriately to obtain the same (3. Sec- 
ond, b is then calculated from a based on equation (|2ip . 
with x depending on i* as given in (jlip . With respect 
to the following discussion we assume here that desorp- 
tion caused deviations of \ from equation (jlip can be 
neglected, which in the framework of the standard nucle- 
ation theory is valid in the long-time limit. In any case, 
more refined estimates of x can easily be implemented. 

As has been shown in equation (|19p a does not depend 
directly on p, but only on (3 and the boundary conditions 
(s). Thus the parameters for the numerical computa- 
tions can be chosen arbitrarily. Using u and t as the 
respective arbitrary units of length and time, the cho- 
sen parameters a all set to unity: s = 1 u, D = 1 u /t 
and F = lpt/(u 2 t). The diffusion field is discretized on 
a 145 x 145 grid for every considered p. To calculate a 
the grid is enlarged by one field in each dimension. To 
evaluate the influence of desorption numerical solutions 
were found for (3 = modeling the no-desorption case, 
(3 = 1.90 which corresponds to the parameters chosen 
for the MC simulation presented in section IIVB 21 and 
finally with (3 = 10 and 50 two cases of medium high and 
very high desorption are discussed. The nucleation prob- 
ability is determined by integration over the stationary 
concentration profiles with exponents i* = 1,2,3,5,10 
and 50 according to equation ©, where 50 is already 
close to the i* — > oo limit and thus the c max ansatz. The 
obtained data is depicted in figure [HI 

Recapitulating the concepts that lead to figure [6] it 
can be read in the following way: if one knows a system 
defined by R-,p, F, D, g, s, A and i* , one can use the cor- 
responding graph to determine by means of b how much 
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FIG. 6: Numerical results for b(s) at different f3 as well as different critical nucleus sizes (i*), which are indicated by the legends. 
The thin dotted lines indicate the expected limit behavior according to equation (|33p for (a) and equation (|35|) for (c) and (d), 
where for (c) only the limit fo(0) is depicted. 



p will change with an infinitesimal change in F at oth- 
erwise constant parameters (especially R). For b > an 
increase in F leads to a longer retained nucleation con- 
trol, while the opposite holds for b < 0. 

All of the graphs in figure [6] exhibit a similar tendency 
of monotonous increase of b with s. The desorptionless 
case is depicted in figure^a). As can be seen, the curves 
for i* = 1,3 and 5 intersect the abscissa meaning the 
direction of the scaling behavior can generally change. 
In the two-dimensional case as s — > the influence of 
the sinks on the concentration field disappears and the 
concentration diverges, which can be seen from the nu- 
merical treatment in the appendix. This does not re- 
quire a divergence of b however. It can also be seen that 
with increasing i* the curves approach the analytically 
approximated behavior for i* — > oo from equation (|33p . 
which is included as the thin dotted line. To yield an ap- 
propriate approximation the corresponding fit parameter 
from table [J in appendix IB1 is used which scales f —>■ Ayr. 
It has to be kept in mind however, that the approxima- 
tion of a radially isotropic concentration field leading to 
equation (|3"3")1 breaks down at large f. 

In figure HJb) one can find the results for (3 = 1.90, 
which correspond to the parameters of the MC simula- 



tions. The obtained data exhibit a significant qualitative 
change in the shape of the curves: at large s the distance 
between the dots becomes smaller which in turn renders 
the 'effective' /3 small; thus the influence of desorption 
becomes less important and the results are very similar 
to the ones obtained in the no-desorption case. However, 
at small s the existence of desorption leads to a distinctly 
different behavior: as before, with s — * the influence of 
the sinks disappears, but now instead of diverging the 
concentration approaches F j g. Consequently the ansatz 
of equation (f34|) holds and the limit for b(s = 0) is given 
by equation (|35|) . In this limit the order of the different i* 
with respect to b is flipped and the curves intersect. The 
final convergence is not shown in the graphs however, as 
it occurs too close to s = 0. 

The approach of the limit at small s becomes more rec- 
ognizable when /3 is raised to 5 as depicted in figure [6]Jc). 
Here it can be clearly seen that with s going to zero the 
obtained data approaches the value of equation (|3"5|) for 
s = 0, which is indicated by the thin dotted lines. Con- 
sidering the other limit [3 is now sufficiently large to show 
an effect even if the effective distance between the dots 
decreases. As this effect also moves the curves toward 
the limit behavior shown by the thin dotted lines in fig- 
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ure [6jd) the intersection is shifted toward 1=1/2. If 
is further increased to 50 (see figure [^d)) the numeri- 
cal data further approach this limit to the point that the 
intersections disappear. The deviations in the region of 
s 0.45 might be attributed to the fact that close to 
coalescence of the dots the area with a low effective dis- 
tance and thus reduced influence of desorption increases. 
For s — > 1/2 the curves seem to approach the theoreti- 
cal prediction again. However, in this region numerical 
artefacts caused by the discretization of the comparably 
small remaining area have to be taken into account, also. 
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2. Simulation 
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Again, to evaluate the theoretical considerations we 
examine the concentration profile around the center of 
the gold dot. To this purpose we simulated two sys- 
tems with n — 20, 30. The obtained radially averaged 
concentrations are given in figure [7] in the dimensionless 
representation. In contrast to the one-dimensional case 
we now consider the radially averaged concentration pro- 
files. The average is taken with respect to only one tile 
of the pattern, meaning from < f < 1/2 the deter- 
mination is straightforward, but for 1/2 < f < 1/V2 
only the corners of the same square contribute. The dot- 
ted lines represent results obtained by numerical solution 
under the given parameters as described in the previous 
section. An approximation formula for this behavior is 
discussed in detail in appendix [5] Of course, the de- 
picted results do not represent a single configuration but 
are the result of an average over the whole simulation. 
This leads to some error with respect to the fact that the 
dot size is not constant, but increases with time via con- 
tinued accretion of new particles. The effective average 
radius chosen for the numerical solution of both problems 
is s c ff = 5.65 er. With p chosen as the geometric mean of 
the substrate dimensions of the MC simulation this leads 
to the reduced radii s c s = 0.27 and 0.18 for n = 20,30, 
respectively. 

As in the one-dimensional case, the predictions based 
on non-interacting adatoms stand in good agreement 
with the simulation results. This is even true for systems 
close to the nucleation threshold: the system with n = 30 
corresponds to the last point of full nucleation control in 
figure[H Also in analogy to the striped pattern, the onset 
of nucleation can be seen from the increase in concentra- 
tion toward the corner of the pattern close to f = 1/ \/2. 
Both cases exhibit a maximal concentration significantly 
lower than the F/g limit, which in dimensionless units 
is given by 2//3 2 = 1.38 and 0.61 for n = 20 and 30, re- 
spectively. These findings correspond to the picture that 
the influence of desorption becomes less important with 
decreasing p and increasing s. 

We further attempt to examine the flux dependence of 
p in the described system. To make use of the previously 
presented numerical analysis (see especially figure [HJb)), 
an estimate of i* is needed. According to the Walton re- 



FIG. 7: Reduced radial concentration for dot patterns with 
lattice spacing n = 20 (o) and 30 (x). The dotted lines rep- 
resent the ideal behavior as calculated by numerical solution. 
Due to the different sizes the reduced dot radii are s e s = 0.27 
and 0.18, respectively. 
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FIG. 8: Nucleation probability in dependence of the reduced 
periodicity for the fluxes 1.43- 10 -7 pt/(step<r 2 ) (o) and 1.43- 
10 -6 pt/(step a 2 ) (x). As before the limit behavior for p <C 1 
and p S> 1 is indicated by dotted lines. The dashed-dotted 
line marks xn — 1/2. The dashed lines symbolize represents 
the estimated expectation in the limit R oc F d according to 
equation (|22a[l using / w 16.5. 



lation all unstable clusters are expected to be in equilib- 
rium while stable clusters continue to grow, thus i* can 
be extracted from the distribution of the cluster sizes. 
This in combination with visual evaluation leads to an 
estimate of i* w 5 for the described system. As can be 
seen from figure IH the regime of nucleation control with 
R <C 1 is found up to p ~ 1.3 which in combination 
with the above A and s c g lead to a reduced dot radius of 
s c ff ~ 0.18. At this s e ff and i* figure [6jb) (which with 
/3 = 1.9 corresponds to the present case) exhibits b w 0.05 
and consequently one expects only a slight F dependence 
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of p. 

With respect to a change in F the Monte Carlo system 
is limited in two ways: on the one hand, a sufficient reduc- 
tion of F requires an according increase in substrate size 
and particle numbers which in turn leads to a significant 
increase in computational effort; while a corresponding 
increase of F on the other hand brings the system to 
its resolution limit, when the observed distances between 
islands reduce to the scale of only a few atoms. 

Based on these considerations a series of simulations 
with a tenfold increased flux of F — 1.43 ■ 10~ 6 pt/(st a 2 ) 
was conducted, the results of which are represented by 
the crosses in figure 0J Similar to before, the three sys- 
tems following the loss of nucleation control were deter- 
mined by twelve independent runs, while the other were 
obtained from six; the nucleus densities were determined 
at the beginning of the saturated regime before the onset 
of coalescence. In contrast to the previous case however, 
this time frame is not as wide as before due to the in- 
creased flux. This leads to a problem, because the scaling 
relations (fTU|) and (fig)) are based on systems at the same 
coverage, not time; hence, for an adequate comparison 
the current simulations would have to be evaluated at 
3 • 10 4 steps, a time at which saturation has not set in, 
yet. An evaluation of both systems at a larger coverage is 
not possible, because of the onset of coalescence. This in- 
compatibility is substantiated by an evaluation of the nu- 
cleus densities. At the increased flux, the surface exhibits 
a higher island density leading to A = 12.9 <r, whereas 
one would ideally expect a value of 10.9 a according to 
equations (fit)]) and (TTTj) with x = 5/7. Additionally, this 
might also be attributed to the comparatively larger size 
of the stable nucleus in relation to A which may promote 
coalescence effects reducing the number of nuclei relative 
to the theoretical expectation and thus resulting in an 
increased value of A. With respect to the loss of nucle- 
ation control the obtained data sets do not exhibit any 
particular difference. A further evaluation based on the 
estimated value of b is not possible, however, as the data 
could not be obtained at the same coverage. 

As an alternative means of comparison between the 
nucleation control of theory and simulation we turn to 
mean number of additional nuclei R, which can be ob- 
tained from via equation ^ and which is depicted in 
figure [5] Its p dependency according to equation (|22a|) 
can be estimated in the following way: taking, as before, 
s w 0.18 for p w 1.3, the numerical evaluation leading to 
figure Ellb) yields a w —3.3. From this one can in turn 
find / R5 16.5 by means of equation (|22b|) with i* = 5. 
Recalling its definition (|22a[) . / is an expression for the 
slope in figure[8] As can be seen from the dashed line, the 
observed behavior corresponds nicely to the theoretical 
estimation, substantiating the validity of the conducted 
ansatz. 



V. DISCUSSION AND SUMMARY 

In this paper we examined the loss of nucleation control 
in epitaxial growth on patterned substrates with special 
emphasis on the pattern spacing. Traditionally the tran- 
sition from (nearly) full to partial nucleation control is 
believed to occur when the nucleus density of the pat- 
tern becomes lower than the density on an unpatterned 
surface (which in a square pattern corresponds to p — 1). 
This is just a first approximation, however. For a more 
detailed understanding we analyzed the formation of the 
first nucleus by means of theory and simulation. 

We conducted continuous Monte Carlo simulations of 
one and two-dimensional patterns and compared these 
results to the theoretical predictions, that are based 
on non-interacting adatoms in the limit of a stationary 
adatom concentration field. In both cases the obtained 
concentration profiles agree well with the predictions, jus- 
tifying the chosen limits which we show are always acces- 
sible at a sufficiently low flux and sufficiently long times. 
Only close to the nucleation of a new dot the mutual 
stabilization of the adatoms becomes relevant. 

Based on the stationary concentration and in principle 
similar to nucleation theor y 13 i 14 i 15 ' 16 we derived an ex- 
pression for the number of additional nuclei (i?) formed 
aside of the pattern by statistical fluctuation. We then 
discuss the scaling behavior (b) of flux (F) and pattern 
spacing (p) at constant R, i.e. xn, by numerical meth- 
ods complemented by analytical evaluations of some ac- 
cessible limit cases. The simulated nucleation behavior 
stands in good agreement with these results: first, as pre- 
dicted , the point where nucleation control is lost does 
not exhibit any significant dependency on the flux; and 
second the expected increase of R with p agrees very well 
with the simulation data. All these observations support 
the validity of the proposed model. 

The present approach, however, does possess cer- 
tain limitations: (a) it is based on the assumption of 
static adatom concentration fields on a patterned sur- 
face, meaning the dot density of the pattern has to be 
similar to the density of a respectively saturated unpat- 
terned surface. We hence require that the experiments 
are quenched at a time where unpatterned surface is 
past the regime of rapid nucleation and in the regime 
of AN /At »0. (b) Due to the assumption of the static 
concentration field the approach is also limited to the 
regime of sn ~ 1- As soon as more than one additional 
nucleus is formed, i.e. xn becomes smaller than 1/2, this 
assumption is not justified any more. This behavior is 
also exhibited by the simulations, (c) Finally, the pre- 
sented approach is not able to explain some behavior 
observed in Ref. [H which exhibited a significant loss of 
nucleation control for p < 1. On the one hand, this suf- 
fered from a lack of experimental data concerning D,g, 
and F. On the other hand the inclusion of dynamic ef- 
fects in a system which has not yet reached the steady 
state might be necessary. 

Therefore, the next step is to generalize the concen- 



12 



tration field approach and to take the formation of more 
than one additional nucleus into account. This can be 
done e.g. by means of a mean field simulation in which 
the discretized concentration field is propagated in time 
and new nuclei are introduced in the sense of a level 
set method. Another important question deals with the 
geometric shape of the clusters. Here comparison with 
free energy minimizations on the basis of density field 
theor y 24 i 25 may yield interesting insight. 
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APPENDIX A: STATIONARY 
CONCENTRATION REGIME 

The following considerations illustrate the limit of the 
stationary concentration field, which provides the basis 
for the calculation of the remaining nucleation probabil- 
ity in section ITVl Starting with the one-dimensional case 
and with g — for simplicity one can express the time 
dependent concentration by means of the limit concen- 
tration 0(^)00 as given by equation (f2"3"| in combination 
with Ac(x,t) describing the difference to c(x)oo: 



c(x, t) = c(x)oo + Ac(x, t) . 



(Al) 



The difference can then be expanded by a Fourier series 
which, having to fulfill the diffusion equation (TT2"1) as well 
as Ac(0, t) = Ac(p, t) = turns out to be 



critical time t c (given e.g. by equation (|A4[) ). This leads 
to a general, monotonously growing dependence of t c (p) 
and its respective inverse p c (t). Second, to provide nearly 
full nucleation control we require a fixed nucleation prob- 
ability in the range < R <C 1, on the basis of which 
the critical flux can be expressed by the general relation 
F c (p,t 1 i?), which decreases monotonously with p and t. 
It should be noted that these critical parameters are not 
to be mistaken for the ones introduced in section [TV] In 
the deterministic limit of i* — > 00 at g = the nucleation 
probability is directly proportional to the maximal con- 
centration from equation (|24[) and one obtains F c oc p~ 2 . 

Finally one can discuss two general scenarios: (a) at 
fixed time scale t the stationary state will be reached for 
all p < p c (t). The second requirement of a fixed R can 
then only be satisfied by a flux of F > F c (p c , t, R). This 
is illustrated by the thick solid line in figure 9(a) (b) on 



the other hand one can consider a scenario at constant 
p. In this case the stationarity requirement is fulfilled if 
t > t c (p) while for R to be retained the flux has to be 
chosen appropriately low satisfying F < F c {p 1 t Cl R); an 
illustration of which can be found in figure 9(b) Con- 



sequently, the requirements of the proposed model can 
in principle be satisfied for any experimental setup if the 
flux is chosen appropriately. 
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FIG. 9: Illustrations for the qualitative examination of F c 
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The Fourier coefficients are determined from the initial 
condition Ac(x, t) = — c(a;) 0o . As the higher modes decay 
with n 2 only the n = 1 mode with 



ai = — tt 
D 



F 4p 2 



(A3) 



dominates the long-time approach of c(x)oo. Concerning 
the exponential decay of Ac(cc, t) one can assume the 
stationary state to be reached when 



t > 



Dir 2 



(A4) 



For the model proposed in section ITVl to hold the con- 
sidered system has to fulfill two requirements: first, to 
reach a steady state the time has to be larger than a 



APPENDIX B: NUMERICAL EVALUATION OF 
THE STATIONARY CONCENTRATION 

In this appendix we evaluate the analytical approx- 
imation of a radially isotropic concentration field from 
equations ([30]) and (|3"Tj) by comparison to the results 
obtained from numerical solution of the diffusion equa- 
tion (fT2|) . The diffusion field is discretized and diffusion is 
expressed by the finite differences of equation (|36p . Sinks 
are modeled by a combination of grid points to represent 
the desired shape. 

First, the hexagonal case, being closer to the circu- 
lar symmetry, is modeled on a 111 x 192 grid with 
Ax = Ay = 1 u approximating the side relation of 
X/Y = 1/V3 at F/D = 40pt/u4 and g = 0/t. A sink 
is placed at the corner and in the center of the unit grid. 
In combination with periodic boundaries this produces 
the desired hexagonal lattice. Three different shapes of 
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FIG. 10: Numerically determined c ma x(s) for circular, squared 
and triangular dots on hexagonal (a) and squared (b) lattices. 
Dot shapes are indicated by the symbols. The dotted lines 
are fitting functions according to equation (|B1|) with the re- 
spective parameters from table [U 



sinks are realized: circles with radius s, squares with edge 
length 2 s and equilateral triangles with edge length 2 s. 

The obtained c max (s) are depicted in figure [TCTa). Ac- 
cording to the theoretical prediction of equation (|3Tj) , the 
data is fitted using 



.fit 



ln(2A s) + 2{A sY 



(Bl) 

with Aq and A\ as fit parameters. Since A\ seems to 
be independent of the dot shape, it is determined from 
an average over fits of the same pattern; Aq is then re- 
calculated based on this result. The obtained fitting pa- 
rameters are given in table |TJ 

The scaling parameter A can be interpreted in the 
sense of an effective radius to approximate a non-circular 
dot by a circle: assuming that for reasonably small s the 
influence of a non-circular dot diminishes rather quickly 



as one moves away from it, 
isotropic profile is obtained; 



pattern: 
dot shape: 



A 
Ai 
B 



hexagonal 
circle square triangle 



an approximately radially 

this could in turn just as 
quadratic 
circle square triangle 



0.938 1.09 0.718 
1.15 1.15 1.15 

73/2 ^3/2 73/2 



0.774 0.913 0.650 
1.28 1.28 1.28 
1/V2 1/V2 1/^2 



TABLE I: Fitting parameters for equations (|B1|) and (|B2|) 
for different patterns and island morphologies. Ao and Ai 
are determined based on numerical solution. B is derived 
geometrically. 



well have been caused by a circular dot of an appropriate 
effective radius. Adopting the picture that primarily the 
area of the dot has to be constant, one can estimate this 
effective radius from simple geometrical considerations. 
For the quadratic and triangular shapes one would expect 
the radii to be related to the circular radius by factors of 
(tt/4) 1 / 2 and (7T/V3) 1 / 2 , respectively. This corresponds 
well to the ratios of the according Aq from table HI 

A square lattice is realized under the same conditions 
as above with the exception of a 135 x 135 grid. The 
adoption of a radially isotropic concentration is not well 
suited for this lattice; nevertheless, the obtained data can 
be properly fitted according to equation (|B1[) as shown in 
figure fTUT b). The parameters were determined as before 
and are given in table fl] Again, the data of all calculated 
shapes can be expressed by the same Ai, but in contrast 
to the hexagonal lattice the ratios of the Aq do not reflect 
the geometrically expected effective radii as closely as 
before. 

Incorporating Aq and A\ into the radial concentra- 
tion (l30l) one obtains 



P so (r,s)=A 1 



{(A Sf - (Br) 



1 



hi 



B f 

Aq 8 



(B2) 

where the introduction of the additional parameter B 
rescales the space coordinate such that the maximum is 
moved from f = 1/2 to the maximal distance of the re- 
spective lattice. B can be derived from simple geometri- 
cal considerations and is also listed in table [H Exemplary 
curves for a quadratic dot of s = 0.07 and a circular one 
with s — 0.15 are depicted in figure [TT] as dashed-dotted 
lines in comparison to the numerically determined pro- 
files. The numerical results (bold lines) were radially av- 
eraged over one unit-cell leading to an edge at f = 1/2, 
where the boundary of the cell is met. This edge is more 
pronounced for the quadratic lattice due to its less cir- 
cular shape. As expected, the maximal concentrations 
(dotted lines) do not exhibit this edge. It can also be 
seen that equation (|B2|) only closely approximates the 
latter part of the profiles as Aq and B do not sufficiently 
compensate to reproduce the used s. 
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FIG. 11: The radial concentration profiles for a circular dot 
with s — 0.07 on a cubic lattice (solid line) and s = 0.15 in 
a hexagonal lattice (dashed line). The bold and dotted lines 
are obtained by numerical simulation and symbolize the mean 
and maximal concentration, respectively. The dashed-dotted 
lines are the approximations according to equation (1B2[) . 
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